GROUND STATES FOR MEAN FIELD MODELS WITH A 
TRANSVERSE COMPONENT 
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Abstract. We investigate global logarithmic asymptotics of ground states for 
a family of quantum mean field models. Our approach is based on a stochastic 
representation and a combination of large deviation and weak KAM techniques. 
The spin- 1 case is worked out in more detail. 



1. The model and the result 

1.1. Introduction. Stochastic representations/path integral approach frequently 
provides a useful intuition and insight into the structure of quantum spin states. 
Numerous examples include |2lf3l[8 , 10 12, 17 18 22l 26] . In this work we rely on 



a path integral approach and related large deviations techniques, and derive global 
logarithmic asymptotics of ground states for a class of quantum mean field models 
in transverse field. These asymtotics limits are identified as weak-KAM 16 type 
solutions of certain Hamilton- Jacobi equations. In principle, such solutions are not 
unique, and an additional refined analysis along the lines of 14 19 20 is needed 
for recovering the correct asymptotic ground state. This issue is addressed in more 
detail for the spin |-case. In particular, our results imply logarithmic asymptotics 
of ground states for models with p-body interactions [5] . 

In the case of Laplacian with periodic potential a weak KAM approach to semi- 
classical asymptotics was already employed in [I]. 

Our stochastic representation gives rise to a family of continuous time Markov 
chains on a simplex (defined below) of The transition rates are enhanced 

by a factor of N, and the chain moves in a potential of the type NF. Ground 
states are Perron-Frobenius eigenfunctions of the corresponding generators. On the 



concluding stages of this work we have learned about the series of papers 23 -25 



The models we consider here essentially fall into a much more general framework 
studied in these works. The authors of 23-25 extend an analysis of Schrodinger 



operators [l4|[l9j[20) on M. d to lattice operators on eZ d , and they develop powerful 
techniques, which go well beyond the scope of our work, and which enable a complete 
asymptotic expansion of low lying eigenvalues and eigenfunctions in neighbourhoods 
of potential wells. 



The paper is organized as follows: The class of models is described in Subsection 1.2 



and the results are formulated in Subsection |1.4| Main steps of our approach are 
explained in Section [2j whereas some of the proofs are relegated to Section [3j The 
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2 DMITRY IOFFE AND ANNA LEVIT 

spin-^ case is studied in Section |4j Finally, in the Appendix, we establish the 
required properties of the Lagrangian C in (1.11) and, accordingly, the required 
regularity properties of local minimizers. 

1.2. Class of Models. Let X be a d- dimensional complex Hilbert space. For the 
rest of the paper we fix an orthonormal basis {|«)} ag _4 of X. We refer to the set A 

of cardinality d as the set of classical labels. Denote projections P a = \a)(a\. The 
induced basis of X^ = CS^X is 

\a) = |oti) <g) • • • <g) | aw) oti, . . . , otjv G A 

The corresponding lifting of the projection operator acting on z-th component is 
P l a = I ® . . . I ® P a ® I g> ■ ■ ■ <8> I. For a G A set = ^Y.i P i- Let M_ N be the 
(^-dimensional vector with operator entries ■ 

We are ready to define the Hamiltonian Hn which acts on X^v, 

-U N = NF (M N ) +J2 B *- (1-1) 

i 

Above, Bi-s are copies of a Hermitian matrix B on X, 5j acts on the z-th component 
of | a). 

We assume: 
Al. F is a real polynomial of finite degree. 

Let Ad be the simplex, A^ = jm G : Y2 m i = m the sequel we shall write 

int (Ad) for the relative interiour of A^. Accordingly, 9A^ = A^ \ int (A^), 

Given m G A^ and a basis vector a G .4.^ let us say that a ~ m, or, equivalently, 
m = m(a), if 

m a = ^ ^ i ^ M^|a) = m a \a), (1.2) 

for all a £ A. Define A^ = A^ PI jf^ d - In other words, m G A^ iff there exists a 
compatible a G A N . In the above notation: 

F (M N ) \a) = F(m(a))\a}. (1.3) 

A2. The transverse field B is stochastic: For any a, (3 G A , 

A Q/3 = A^ ^ (a|B|/9) > 0. (1.4) 

Furthermore, A is an irreducible kernel on A. Without loss of generality we shall 
assume that A = on the diagonal. 

1.3. An Example: Spin-s Models. The relation between the dimension d of X 
and the half- integer spin s is d — 2s + 1. The set of classical labels is 

A = {— s, -s + 1, . . . , s} . 

The stochastic operators are B{ = ASf . A > is the strength of the transverse field. 
Altogether, the Hamiltonian is 

-U N = NF (M5 S , M5 s+1 , . . . , M S N ) + A J2 S *- (1-5) 
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For instance, the case of p-body ferromagnetic interaction corresponds to 

F(M»,M»_ V ...,M»)= fe^X = (±^S'X . ( L6 ) 
The operators S x act (under convention that |s + 1) = | — s — 1) = 0) on X as 

S» = ^(s + 1) - «(« - 1)|« - 1) + -^(s + 1) - a(a + l)|a + 1) (1.7) 

_ — 

Consequently, the jump rates X a s are given by 



v/s(s + 1) - a/3, if|a-/3| = l 



0, otherwise. 



1.4. The Result. In order to develop an asymptotic description of finite volume 
ground states we need to introduce some additional notation: For m G A^ set 



c N (m) 

The vectors \m) G Xn, 



N \ N\ 



NmJ H{Nm c 



) = ~r=^= E la> 



are normalized and orthogonal for different m G A 



N 



By Perron- Frobenius theorem and Lemma 2.1 below the ground state of Hn is 
fully symmetrized, that is of the form 

\h N ) = h N(m)\m), (1.9) 

meA^ 



and h^{rn) > for every m G A^ (see Subsection 2.1). Let us represent 



h N (m) = e ~ N ^^\ (1.10) 

It would be convenient to identify ip^ with its linear interpolation (which is an 
element of the space of continuous functions C(A d )). 
Next introduce: 

"Ho {m, = E \ /ma7n f } (cosh(9f} - 9 a ) - 1) 

Q ^ (1.11) 

£o (zzi, u) = sup { (v, 8) — Ho (m, 9) } 
e 

For m G Ad define 

^(m) = ^X) A ^(V^- V^) 2 -^)- (1-12) 
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Finally set X a = J2p Kp, 

U (m, 6) =H (m,6) - V(m) 

= ^2 y/m a m^X a pco8h(6fi - 9 a ) - ^m Q A Q + F(m) 

a/3 a 

and 

C(m,v) = Co(rn,v) + V(m). 
In (1.13) above (•, •) is the usual scalar product on M d . 



;i.i3) 



Theorem A. Let E\j be the bottom eigenvalue ofH N . Set X = ^2 a X a . Then the 
limit 

-A + n= hm ^ (1.14) 



exists. Moreover, 



iV-s-oo N 

X\ = min V(m). 



;i.i5) 



Furthermore, the sequence {iPn} is precompact in C(A d ). Any subsequential limit ip 
satisfies: For any T > and any m e A^, 



?p(m) 



inf L( 7 (0))+/ £(#),7 / W)dt-Tr 1 ), (1.16) 

7 : 7 (T)=m ^ ,/ J 

where the infimum above is over all absolutely continuous curves 7 : [0, T] — > Ad- 

Moreover, the set Ai^ of all local minima of ip is a subset of Aiy — axgmin(V) C 
int(A d ). 

Remark 1. Hamiltonians Ho,!! are invariant under the shifts 9_ h-> 6 + cl, and, as 
a result, the Lagrangians Cq,C are infinite whenever (y_,l) 7^ 0. Also, Co(m,0) = 
= mm£ (m,v). Consequently, C(m,0) = V(m) = — %(m, 0), and (1.15) could be 
rewritten as 

X\ = min C(m, v) = — max'H(m, 0). (1-17) 



Remark 2. Either of (1.15) and (1.16) unambiguously characterizes \\, but not ijj. 
As we shall explain in the sequel, if ip satisfies (1.16), then the weak KAM theory 
of Fathi fld^ implies that ip is a viscosity solution (see Subsection 2.1 for the precise 
statement) on int (Ad) of the Hamilton-Jacobi equation 

H (m, Vip(rn)) 



1.18 



Note that since if) is a function on Ad, the gradients Vip lie in the subspace 

R d o = {v:(v,l) = 0}. (1.19) 



In general there might be many viscosity solutions of (1.18) which comply with the 
conclusions of Theorem [A] The solutions which are subsequential limits of {^n} are 
called admissible. Although we expect uniqueness of global admissible solutions for 
a large class of models, our approach does not offer a procedure for selecting the 
latter. The viscosity setup is important - at least for a large class of symmetric 
potentials the global admissible solutions are not smooth and develop shocks. A 
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proper selection procedure should be related to a more refined analysis of the low 
lying spectra of Hn, As it was mentioned in the Introduction sharp asymptotics 
of eigenvalues and eigenfunctions in vicinity of potential wells were derived in a 
much more general context in [23TI25]. In particular, it is explained therein how 



such asymptotics are related to (smooth) local solutions of (1.18). Implications of 
these results for a characterization of global admissible solutions is beyond the scope 
of this work and hopefully shall be addressed in full generality elsewhere. In the 
concluding Section [i] we work out a particular case of spin-| models.. 

2. Structure of the theory 

2.1. Spectrum of Hn- Let H S N be the sub-space of X^ which consist of those 
vectors \b) which do not vanish under symmetrization. Namely, \b) = Y2 a a a\^) 
X S N \ 0, if a nQL 7^ for some a G A N , where 7r is a permutation of {1, . . . , N} 
with ira(i) = afc). The sub-space X^ is invariant for the Hamiltonian Hn- The 
ground state of always belongs to ~K S N . For the rest of the paper we shall work 
with the restriction of T-Ln to ~K S N . 

All eigenfunctions of Hn (restricted to X^) have mean-field representatives: 

Lemma 2.1. If E^ is an eigenvalue ofH^, then there exists a function on 
such that \h N ) — J2 m eA N ^n(ir) \WlI ^ s a corresponding eigenfunction: 

U N \h N ) =E N \h N ). (2.1) 

Proof. Let En be an eigenvalue of %n- Let |&jv) = J2 a &A N a oJ^) e ^-n be an 
eigenfunction corresponding to the eigenvalue En- Let C(a, 0) = ((3\B\a) be the 
matrix elements of B = Yli^i- Thus, B\a) = C(a, The eigenfunction 
equation is recorded as: V/3 

22 aqC(a, 0) = {-E N - F(m)) a p. 

a 

Note that C(a,/3) = C{^a^0). Consequently, since in addition m{0) = m(n(3), 

22 a ™C(a, 0) = (-E N - F(m)) a n p. 

a 

Therefore, \kon) — J2 a a ^\^) * s a ^ so an eigenfunction. Since the sum X^ ™ does 
not change if we permute the entries of a, and since, by assumption \bx) G X° N , the 
claim follows with \hpf) — l 71 "^)- D 

2.2. Stochastic Representation. Let a(t) be the continuous time Markov chain 
on A with jump rates A Qj g. is the path measure for N independent copies of 
such chain starting from a. Then the following representation of the entries of the 
density matrix holds [3 22: For any T > and any a, G A N 

e- NXT (0\e- Tn -\a) = Ef exp jiV £ F(m(t))dt| l {a(T)= ^ } . (2.2) 



DMITRY IOFFE AND ANNA LEVIT 



Above m(t) = m(a(t)). 



2.3. Mean Field Lumping. The process m N (t) = m(t) = is a continuous 

time Markov chain on A^ with the generator 

QNf{m) = N ^ m a \ a p If I 
It is reversible with respect to the measure 



m 



8p - 5 a 
N 



f(rn) 



(2-3) 



A c N (m) 



(2.4) 



Summing up in (2.2), 



- NXr (m'\e- TH 



\m) 



fi N (m 



'^E^exp^N^ F(m(t))dtJ l {m{T )= m >}, (2.5) 



for every T > and every m,m' E A% . 

Using Girsanov's formula one can rewrite (2.5) in a variety of ways for different 
modifications of the jump rates in (2.3). Namely, let g be a positive function on 
A^j . Consider the modified rates 

Let P^' 9 be the corresponding path measure. Then, the right hand side of (2.5) 
reads as t 



where 

F g (m) = ^ Kf) m c 

a/3 

A self-suggesting choice is 
g(m) = 



g{m + {5 p -5 a )/N) i 
g(m) 



{m(T)=m'} ; 



+ F(m). 



(2-6) 
(2-7) 



'fjLN{m) 



\ N J = N^m a {m p + l/N) \ a/3 . (2.8) 



For the rest of the paper we fix g as in (|2.8|). The corresponding generator 

5p - 5 a 



N 



f(m] 



(2-9) 



is reversible with respect to the uniform measure on A^. The function F g in (2.7) 
equals to 

F g (m) = -V(m) + Z N (m), (2.10) 
where V is precisely the function defined in (1.12), and the correction 



a/3 



(2.H) 
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All together, (2.5) reads as 
e-^^'le-^I^^E^expj-AT^ (V - E N ) (m(t))dtj t {m{T)=m , } , (2.12) 

As we shall see below it happens to be convenient to work simultaneously with both 
representations (2.5) and (2.12). 



Note that an immediate consequence of (2.12) is: 

Lemma 2.2. En is an eigenvalue of Hn with [itjv) — Y1 u n('D1)\ir) being the cor- 
responding normalized eigenfunction if and only if un is also an eigenfunction of 



s N = g 9 N + NF, 



G 9 



N(V - E N ) with 



R 



N 



(NX + E N ) 



(2.13) 



being the corresponding eigenvalue. 



2.4. The Eigenfunction Equation. Assumption A. 2 and Perron-Frobenius the- 
orem imply that %n has a non- degenerate ground state |/ijv) = Ylim ^Jv(2Zl) Izzi) with 
strictly positive entries /iAr(m) > 0. By Lemma 2.2, hjsr{m) is the principal eigenfunc- 
tion of G% + NF g (m) with the corresponding top eigenvalue —R N = — (NX + E N ). 
The corresponding eigenfunction equation is: For any T > 0, 



h N (m) 



E^expj-A^ (V-E N )(m(t))dt + TR N ^h N (m(T)), (2.14) 



By reversibility, 

Lemma 2.3. Functions {h^} satisfy: For every T > 0, N and m e A^ 



h N (m) = h N (rnf)E^ex P {"^ (V - E N ) (m(t))dt + TR 1 ^ 1 



{m(T)=m}- 

(2.15) 



2.5. Compactness and Large Deviations. For m,m' e A^ define 



Z^ ,9 (mf, m) 



logE^exp { -N I (V- ~ N )(m(t))dt } 1 



N 



I{m(T)=m}; 



(2.16) 



In the sequel we shall identify Z. 



N,g, 



with its continuous interpolation on A^ x A^. 



Let ACt be the family of all absolutely continuous trajectories 7 : [0, T] i-> A<f. 
For m, m' G A^ define 



ZtL(m', m) 



inf 

l(0)=m'MT)=mJ 
■y&AC T 



£( 7 (t),7'W)dt, 



(2.17) 



where the Lagrangian C was defined in (1.13). 
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rN,g 



IS 



Theorem 2.1. For all T sufficiently large the sequence of functions |Z T 

equi- continuous on A d x A^ and uniformly locally Lipschitz on int (Ad x A^). Fur- 
thermore, for all T sufficiently large, 



lim 

Af->oo 







(2.18) 



Z^ ,9 (mf ,m) — Z!f(m,m) 
simultaneously for all m,m' G A^. 

Note that the equi-continuity of in Theorem 

gence in (2.18) is actually uniform. Consequently, Z^(-, •) is continuous on A^ x A 
and locally Lipschitz on int (A d x A^). 

Theorem 2.1 is a somewhat standard statement. Its proof will be sketched in 



2.1 



implies that the conver- 



Subsection |3~T 



2.6. Lax-Oleinik Semigroup and Weak KAM. Recall the representation of the 
leading eigenfunction h^(m) = e~ N ^ N ^—\ In the sequel we shall identify ipN with 
its (continuous ) interpolation on A d ; ip N G C(A d ). 

Theorem 2.2. The sequence of numbers R l N /N is bounded in 1R. The sequence of 
functions {4>n} is equi- continuous on A d and uniformly locally Lipschitz on int (A d ). 

Proof. Since is the Perron-Frobenius eigenvalue, 

"^T^^ logE 5 9eXp { N J F ^( s )) ds }' 



R N 

N 



which is bounded since F g is bounded on A<j. 



On the other hand, by (2.15), the 



equi-continuity and the uniform local Lipschitz property of {iPn} is inherited from 



the corresponding properties of < Z, 



□ 



Proof of (1.14) and (1.16) of Theorem [Aj Theorems 2.1 2.2 and Lemma 2.3 
imply that any accumulation point (r, ifi) G 1R x C(A d ) of the sequence {-^-Rjy, i/jn\ 
satisfies: ip is locally Lipschitz on int (A^) and 



ib(m) = inf 

7 (T)=m 



^ (7(0)) + 



J £(7(t),7'(t))dt| -Tr = U T ^(m) - Tr. (2.19) 



for every T > and each m G A^. Accumulation points of ip^ are called admis- 
sible solutions of (2.19). Since Ut is non-expanding on C(A^), validity of equation 



(2.19) unambiguously determines r, which implies that the limit ri = lim 



R N 

N—too ]\f 



indeed exists. In view of (2.13) we, therefore, have established (1.14) and (1.16) of 
Theorem El □ 
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2.7. Viscosity Solutions. Recall the definition (1.13) of T-L(m,6). For r x defined 
as above consider the Hamilton- Jacobi equation 

H(m, Vip(m)) = -r x . (2.20) 

Definition 1. For m G A^ and ip G C(A^) define lower and upper sub-differentials 



D_ip(m) 



lim inf 4>(™!)-4>(m) - (£,m'-m) > Q 



\m — ml 



(2.21) 



and, similarly, D + ip(m) with lim inf changed to lim sup and the sign of the inequality 
flipped. 

A locally Lip schitz function ip is said to be a viscosity supersolution of (2.20) at 
m ifH(rn,£) > —X\ for every £ G D^ip(m). Similarly, it is said to be a viscosity 
subsolution of \2.20 ) at m ifH(m,£) < — ri for every £ G D + ip(m). ip is viscosity 
solution of \2. 20 ) at m if it is both sub and super viscosity solution. 



Theorem 2.3. If ip is a weak-KAM solution (of (2.19) with r = r ly ) ; then it is a 
viscosity solution of (2.20) on int (A rf ). 



The proof of Theorem |2.3| is relegated to Subsection |3.2 
2.8. Minima of ip. 



Theorem 2.4. Let ip be a weak KAM solution of (2.19). Then all local minima of 
if) lie in the interiour int (A^). 



Theorem 2.4 will be proved in Subsection 3.3 



2.9. Stochastic Representation of the Ground State. The eigenfunction equa- 
tion (2.14) defines a Markovian semi-group 



E$f(m) 



h N (m)~ 



E% 9 exp \ N / F g (m(t))dt + TR l N \ h N {m{T))f{m{T)) 



This corresponds to continuous time Markov chain with the generator 



(2.22) 



G 9 N f{m) 



h N (m) 



a,/3 



5p - So 
N 



f lm + 



$1 - So 
N 



fitR, 



(2.23) 



In the sequel we shall refer to Q g N as to the generator of the ground state chain. 

Lemma 2.4. The generator Q 9 N is reversible with respect to the probability measure 
u^r(m) — h 2 N (m) = e~ 2N ^ N ^ m \ Furthermore, En is an eigenvalue ofH-N if and only 
if Ejf — En is an eigenvalue of Q 9 N . 

It is straightforward to check that Q 9 N satisfies the detailed balance condition 
with respect to u N . It is equally straightforward to see from (2.22) that g^ is an 
eigenfunction of Q 9 N + NF g , and hence by Lemma 2.2 of Hn, if and only if gN/hN 
is an eigenfunction of Q 9 



N ' 
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3. Proofs 



3.1. Proof of Theorem 2.1 Consider the family of processes {m(-) = m Ar (-)} with 
generator Q 9 N defined in (2.9). We shall identify m N with its linear interpolation. 
For each T > 0, the family {m N (-)} is exponentially tight on Co,t (A<z)- 

Recall the definition of Co in Q1.11J ). One can follow the approach of 15 and 
to combine the Large Deviation Principle for projective limits fl3] with the inverse 



contraction principle of [21] in order to conclude: 

Lemma 3.1. For each T > and every initial condition m G A^ (where for each N 
we identify m with its discretization \ Nm\/N G A^ ) the family of processes {m(£)} 
satisfy a large deviations principle on C^t ( A^) with the following good rate function 

J Q T C ("y(s), 7'(s))ds, if ^ is absolutely continuous 
J T (7) = { and 7 (0) = m . (3.1) 

otherwise. 

By the upper bound in Varadhan's lemma, 

limsup Z^' 9 (mf, m) < Z^,(rn , m) 

On the other hand, by the lower bound in Varadhan's lemma, for each 5 > 
Z^(m',m) < liminf sup Z T ' 9 (m', m-J. 



Therefore, (2.18) is a consequence of the claimed continuity of m — > Z^' 9 (rnf ,m). 
Let us proceed with establishing the asserted continuity properties of the family 
Z^' 9 (-, •). By reversibility, 

Z^' 9 (m,m) = Z^ ,g (m' , m), (3-2) 

so it would be enough to explore those in the second variable only. 
An equivalent task is to check continuity properties of 

Z^{mf,m) = llogE^exp jiV^ F(m(t))d* j l {m{T )= m }, (3-3) 

Indeed, by (ggb 



Zr' 9 (m, m') = \ / — rZSim, m!) 

K J y Mm') Tl -'- ; 

Now, under the process m N {t) is a super-position of A" independent particles 
which hop on the finite set A with irreducible rates \ a p. Since F is bounded, the 
following claim is straightforward: 



Lemma 3.2. There exist To > 0, 6q > and a constant c\ < oo such that 

e NZ^_ e (rn',m) > ^Ne^NZ* (m' ,m) _ /g ^x 

uniformly in N , T > T , e < e , rnf,m G A rf and A". 
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Indeed, trajectories m(-) on [0,T] and trajectories m(-) on [0, T — e]: are related 
by the following one to one map: fh{t) = m(t^-). Since for some c 2 = c 2 (T ) > 0, 
up to exponentially small factors, the total number of jumps of all the particles is at 



most C2NT, the Radon - Nikodym derivative is under control and (3.4) follows. □ 
As a result, for any T, e as above, and for any mf, m 1 , m 2 G A^, 

Z^(wf,m 2 ) > Z^(mf,m 1 )-c 1 e + Zf(m},m 2 ). (3.5) 

Fix ra^ra 2 and define A+ = A+(rn},rn?) = {a : m\ > m 2 }. For a G A+ define 
5 a = m\ — m 2 . One way to drive m(-) from m 1 to m 2 during e units of time is to 
choose NS a particles out of Nm} a for each a G A + , and to redistribute them during 
e units of time into A \ A + , without touching the rest of the particles. There is an 
obvious uniform lower bound c^e n that a particle starting at the state a will be at 
state p at time e. We infer: 



OS.4+ 

Hence, 

/7JV 



c (m\m 2 ) > -c 4 e - c 5 ^ 5 Q fdlog * - log ^ . (3.7) 
Both, the equi-continuity of m — )• i^Y (m', m) on A^ and its uniform local Lipschitz 



property on int (A^) readily follow from (3.5) and (3.7). □ 



3.2. Proof of Theorem 2.3 We follow the approach of |16|: Let m G int(A) 
and assume that u is a smooth function such that {m} = argmin {u — ip} in a 
neighbourhood of m. Then, 

u{m) < M ( 7 (-t))+ / £(7, 7 ')ds - nt 



for any t > and for any smooth curve 7 with 7(0) = m. Let v = 7'(0). Then, 

Vtt(m) • v — C(m,v) < — r%. 

Since the above holds for any m G Rq, 'H(m, Vu(m)) < — ri follows. 

In order to check that ^ is a super-solution, note that by the upper and lower 
bounds on the Lagrangian £ derived in the Appendix, and by the local Lipschitz 
property of (bounded and continuous) ip the minimum 

V>(7(0))+ / V(7,7')-ri)d S 



mm 

7(to)=m [ ' " ' " J Q 

is attained at some 7* with 7 *(0) = m' in a (^-neighbourhood of m, for all to and Sq 
appropriately small. As it is explained in the Appendix, the minimizing curve 7* is 
C°° and stays inside int (A^). Evidently, 

ip(m) = V(7*(*)) 



y ° (£(7*, 7i) - rx) ds 
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for every t G [0, to]- Assume that u is smooth and argmax{w — ip} = {m} in a Sq 
neighbourhood of m. Then, 



u (m) = u(7*(t )) > u Mt)) 



for every t G [0, to]- Set v = 7* (to)- We infer: 

V«(m) • v > C(m, v) — X\. 
Consequently, "H(m, V«(m)) > — ri. 



□ 



3.3. Proof of (1.15) of Theorem A and Theorem 2.4 By Theorem 2.2 and 
since v N {m) = e r 2JVl M^ ) it follows that 

nnn^(m) = 0. 



Let us rewrite (1.16) as 



^(m)= inf V(T(0))+ / (£(7,7')-ri)dt 

7(T)=m ^ Jo 



(3.* 



Since V(m) = £(m, 0) < £(m, v), the above might be possible only if x x = min^ V(m). 



Furthermore, the Lagrangian £ is uniformly super-linear in the second variable: 



By (A.l ) of the Appendix for every C > and 5 > we can find T > such that 

inf / T (£(7,7')-ri)dt>a 

diam(7)>5 



Which means that for C > maxip, the contribution to (3.8) for 7-s with the diameter 
larger than S could be ignored. 



Let, therefore, diam(7) < 5 . By (1.17), 

/ (£(7, 7') - rx) dt > T min (V(mf) - 
Jo rnl&i 



We infer: 



•0(m) > min ip{w!.) + ^ mm (Viwl) ~ r i) 

\m'—m\<S \m'—m\<5 



(3.9) 



(3.10) 



The claim of Theorem 12.41 follows as soon as we notice that all the minima of 
m 1 — y V(m) belong to int (A^). □ 



4. Results for Spin-^ Model 



For spin-s models (1.5) 



i"i = min V(m) = min < — v/ s(s + 1) — a(3 {Jvn a — Jvtlr) 2 — F (m) 



\*-P\=l 



(4.1) 
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In spin-| case it is convenient to take { — 1,1} instead of { — |, ^} as a set of 
classical labels for Spin-^ Model. The Hamiltonian is given by 



-n N = nf (m^i, m n ) + \J2 °i 

i 

where 

*-=(?j) ^(j^). (4.2) 

In this notation a z \oc) = at\ct) and a x \a) = \ — at) for a = ±1. 

The simplex is just a segment {(m_i, mi) : m_i + mi = 1}, parametrized by a 
single variable m = mi— m_i G [—1, 1]. Any vector 9 G M 2 , * s °f the f° rm £ — —9) . 
Define F{m) = F (^-^, Thus, in terms of m and 0, the Hamiltonian in (1.13) 

is 

■H(m, 9) = Wl - m 2 cosh(2#) - A + F(m). (4.3) 
Consequently, the effective potential 

V(m) = -H(m, 0) = A - (\y/l - m 2 + F(m)) , (4.4) 

and the asymptotic leading eigenvalue ri is given by 

n = min V(m) = A — max < A\/l — m 2 + F(m) > . (4.5) 

me(-l,l) me[-l,l] I J 



4.1. Minima of V . In order to explore the minimization problem (4.5) it would be 
convenient to represent F{m) = —G (sign(m)\/l — m 2 ). Then, 

ri = A - max {X\t\ - G(t)} . (4.6) 

-KKl 



Indeed, as it clearly seen from (4.5) (and as it follows in general by Theorem 2.4), 
all the minima of V belong to (—1, 1), a possible jump discontinuity of G at zero 
(if F(—l) ^ F(l)) plays no role for the computation of maxima. Note also that 
G(-1) = G(1) = F(0). 
Let 

M\ = argmim (V) = {m G (—1, 1) : %{m, 0) = — r x } . 

In other words, let T\ C [—1,1] be the set of maximizers in (4.6). We set T\ = 
U T A ~, where T A + = T A n (0, 1). Then, 

m G M x t = sign(m)Vl — m 2 G 7^ ign(m) . (4.7) 

By assumption Al, G is a polynomial of finite degree on each of the intervals [—1, 0) 
and (0, 1] , so the set 7a is finite, and its maximal cardinality is deg(G) — 1. Various 
options are depicted on Figure [l] (for simplicity we depict only the (0, 1] interval 
and, accordingly, the set T x + ): 

(a) First of all there exists A c G [0, oo), such that 7a = {±1} on (A c , oo). 

(b) If A c > 0, then 7a c still contains ±1. It could happen, however, that 7a c = 
{— l,ii, . . . ,tk, 1} contains at most k < deg(G)/2 other points. In the latter case 
(—1, 0) U (0, 1) necessarily contains at least 2k inflection points of G. 
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G[t) 




G{t) 



Figure 1. (a) The critical A c and = {t\, 1}. (b) T x + is a single- 
ton. 

(c) T x + = {^1,^2,^3}- There are at least two inflection points of G on 
(*i,ta) 

(c) There might be other exceptional values of A < A c for which either of is not 
a singleton. If, for instance T x + = {t^, . . . , i£} is not a singleton, then the interval 
(ti,tk) contains at least 2(k — 1) inflection points of G. Since there are at most 
deg(G) — 2 inflection points all together, and since intervals spanned by different 
7^ are disjoint, we infer that is not a singleton for at most deg(G)/2 values of 
A. 

Values of A for which the cardinality of 7a changes correspond to first order phase 
transitions in the ground state. 

4.2. Ferromagnetic p-body interaction. In the usual Curie- Weiss case with pair 



interactions G(t) = | (t 2 — 1), so that ri 



(A-l) 2 



if A < 1, and, accordingly, r x = 



if A > 1. For A > 1 the set M x = {0}. For A G (0, 1), M x = {±Vl - A 2 }. No first 
order transition occurs. 

In the p > 2-body ferromagnetic interaction case ( |1.6[ ) the function 

G(t) = -sign(ty{l-t 2 ) p/2 . 

For odd p maximizers of X\t\ — G(t) always lie in (0, 1]. For even p the set Ai\ is 
symmetric. Thus in either case it is enough to consider 

r^A-maxjAt+fl-ty 72 ). (4.8) 
te(o,i] I v ' J 

The crucial difference between the Curie- Weiss case p = 2 and p > 2 is that in the 

latter situation, G'(l) = 0, and G contains an inflection point t p 

(0, 1). An easy computation reveals that for p > 2, 

1 



— l -r inside 
p-i 



A, 



V 



p — 1 

Accordingly, for even p, 



(p-1) 



£-1 
2 ± 



and Tx 



p — 1 



M Xc = {0,±m} = {0,±< 



I ' pip -2) 
(p-l) s 



(4.9) 



(4.10) 
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whereas for odd p; M.\ c = {0, m} = jo, y | ■ This is precisely formula (14) 
of |5). Fo r A > A c the set _M A = {0}. For A < A c there exists m* = m*(X,p) G 

(yjp^ip'' l) sucn ^ ne se ^ * s a s i n gl e t° n {m*} in the odd case, whereas 
Aix = {±m*} in the even case. Thus, for mean-filed models with p-body interaction, 
A c is the only value at which fist order transition in the ground state occurs. 

4.3. Asymptotic ground states. Let us return to general polynomial interactions 
F. Fix m G (—1, 1) and consider the equation, 

H{m,&) = -r l . (4.11) 

The Hamiltonian %{m, •) is strictly convex and symmetric. Hence, if m G then 
6 = is the unique solution. If m G" A4\, then necessarily %(m,0) < — rj.. Hence, 
there exist 9{m) > 0, such that ±9(m) are the unique solutions to (4.11). If F is 
symmetric, then 6(—m) = —6{m). In any case, however, the following holds: 

Let ip be an admissible solution of "H(m, ip') = — K\. Since ip is locally Lipschitz, it 
is a.e. differentiable. Consequently, ifj'(m) — ±0(m) a.e. on (—1,1). The proposition 
below relies only on the fact that ip is a viscosity solution on (—1,1). 

Proposition 4.1. Record Ai\ = { — 1 < . . . , < 1} in the increasing order. 
Set m = —1 and = 1. T/ien on eac/i of the intervals [m^, m^ + i] £/ie gradient 
ip' is of the following form: There exists m\ G [m^, m^+i], stzc/i that: 

ip' = 6 on [m^,m|) and ■0' = — 6 on (m|,m^ + i]. (4-12) 

Proof. It would be enough to prove the following: If m G (mt,m£+i) and ip'(m) = 
0(m), then for any n G (rri£,m), 

rm 

ip(n) = ip(m) - 9(t)dt. (4.13) 
Recall that since z/> is a viscosity solution on (—1,1), then 

liminf + 7 ~ >e^9£ (-e(n*),e(n*)). (4.14) 



for any n* G (—1, 1). We shall show that if (4.13) is violated for some n G (m£,m), 
then (4.14) is violated as well in the sense that there exists n* G (n, m) and # G 
(—0(n*),6(n*)) such that the right hand side of ( |4.14[ ) holds. 
Indeed, since ^'(O = a - e - on ( — 1? 1) it always holds that 



V>(n) > if)(m) - / 0(t)dt 

J n 

Let us assume strict inequality. For k G [n, m] define 

By construction v(n) = 0, t> (m) = p(m — n) < (m — n) and v'(m) = 1. There is no 
loss of generality to assume that p > 0. Hence there exists n* G (n, m) such that 
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v(m) 




Figure 2. p e dv(n*). 

p E dv{n*) (see Figure [2]). By continuity of 9{m) this would mean that 
^(n* + e)-V(n*) 



lim inf 



> p6(n*) - (1 - p)0(ra*) = (2p - l)9(n*). 



Since for any inner point n* G (m^,m^ + i); #(?2*) > 0, we arrived to a contradiction. 

□ 



Remark 3. Note that Proposition 4-1 implies that ground states if) with more than 



one local minimum necessarily develop shocks. 



2.4 



in the latter case, if)' = —9 on 



If me < m\ < m^ + i from Proposition 4.1[ then m\ is a local maximum of if), and 
m\ is a shock location for the stationary Hamilton- Jacobi equation T-L{m,ip) = —X\. 

If ip is, in addition, a weak KAM solution (in particular, if ip is admissible), then 
argmin {ip} C Consequently, by Theorem 

(—1, mi) and ip' = 9 on (m^, 1). 

Admissible solutions are always normalized in the sense that minip = and the 
minimum is attained on Ai\ C (—1,1). It follows that admissible are uniquely 
defined in the following two cases: 

CASE 1. The set M\ = {m*} is a singleton. Then, if/ = —9 on (— l,m*) and 
if)' — 9 on (m*, 1). Consequently, 



if){m) 



9{t)dt 



(4.15) 



CASE 2. The interaction F is symmetric and Ai\ = {±m*}. Then, if) is also 
symmetric; if)' = 9 on (— m*, 0) U (m*, 1) and if)' = —9 on (—1, 0) U (0, m*). That is 
ij)(m) = ip(—m) and if) is still given by (4.15) for m > 0. Note that in this case if)' 
has a jump at m = 0. 

Let A c be the set of A which do not fall into one of the two cases above. As we 
have seen in Subsection 4.2, A c = in the case of Curie- Weiss model, and A c = {A c } 



see (4.9)) for general p-body interaction. 



STOCHASTIC REPRESENTATION OF GROUND STATES 



17 



4.4. Multiple Wells. We shall refer to A G A c as to the case of multiple wells. 
Note first of all that there is a continuum of normalized solution of (1.16) as soon 
as the cardinality > 2. Indeed, it is easy to see that any normalized ip which 

complies with the conclusion of Proposition 4.1 will be a solution to (1.16). 

One needs, therefore, an additional criterion to determine locations of shocks {m^} 
or, equivalently, to determine values {ijj(me)} for admissible solutions. It would 
be tempting to derive location of shocks by some natural limiting procedure via 
stabilization of shock propagation along Rankine-Hugoniot curves. Since however, 
we arrived to ( 1.16 ) directly from the eigenvalue equation without recourse to a finite 
horizon problem, it was not clear to us which limit to consider. Our selection of 
admissible solutions to (1.16) is based on a refined asymptotic analysis of Dirichlet 
eigenvalues in a vicinity of points belonging to the set M.\ . Namely, a point 
me G M.\ can be local minima of an admissible solution if> only if there is an 
exponential splitting of the corresponding bottom eigenvalues. Precise result is 
formulated in Proposition 4.2| below. 

enable to explore asymptotic expansions of such eigen- 



The results of 23, 24 , 24 



values with any degree of precision. In the simplest case we deduce the following: 
Corollary 4.1. Assume that 

A 



mm Xo{ m ) 

m£M\ 



mm 



1 



77T 



Vl - m 2 F"(m) 



(4.16) 



is attained at either a unique point m* (non- symmetric potentials) or at a unique 
couple ±m* (symmetric potentials). Then there is a unique admissible solution ip, 
which is still given by ( |^.15| ). 

For instance, in the critical (A = A c ) case of p > 2 body interaction, a substitution 
of {D) and flUu) yields: 

Xo(0) = A c and 



Xo(m) = (p- 2)(p - 1)A C . 



(4.17) 



In other words, Xo(0) < Xo(w), for any p > 2 and A = X c (p) Consequently, even 
at A = A c there is still a unique admissible solution ip(m) = \f™ 9(t)dt\ with the 
unique minimum at m* = 0. 

We explain Corollary |4.1| in the concluding paragraph of this Section. 



Spectral Asymptotics and the Set M.\. 

denote M. = {mi, . . . m^}. 



Assume that A G A r and, as before, 



Lemma 4.1. For any 5 > there exists e > such that 

ipijn) > e, 



mm 

d(m,M\)>S 



(4.18) 



uniformly in normalized admissible solutions of (1.16). 
Proof. Let m G {m h mi + i) 

ijj(m) > min 



By Proposition 4.1 

i/j(mi) + 



0{t)dt , ^(m w ) + 



e(t)dt 



and (4.18) follows. 



□ 
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; S^h 



N 



DMITRY IOFFE AND ANNA LEVIT 

-Ni/> N 

—R^tiN- Recall: 

S N f(m) = N(F(m)-\)f(m) 

NX 



+ 



+ 



2 
NX 



(1 - m)(l + m + — )/ I m + 
(1 + m)(l - m+ -)/lm- 



AT 
2 

N 



(4.19) 



Pick < 5 < ^min/|m; + i — mi\. Let 1 = YloXi De a partition of unity satisfying: 



For I = 1 jfe 



Xi = 1 on 7 5 (mj) and Xi = on I 26( m i)- 



Above I^m) is the interval [m — rj,m + rj\. By Lemma 4.1 there exists e > such 
that for I = 1, . . . , k and all A^ large enough 

^ log max| (S N + Rjf) Xih N {m) \ < -e. (4.20) 

Let 5^ be a Dirichlet restriction of S to Is (mi). Let —R^i be the leading eigenvalue 
of 

We are entitled to conclude: There exists e' > such that 

1 



N 



log\R N — mini^jy;! < — e'. 



(4.21) 



Furthermore, 



Proposition 4.2. If I = 0, . . . , k and ip = lim^oo ip N . is a subsequential limit such 
that mi is a local minimum ofip, then there exists e' > such that: 

^log|i4.-i4.J<-e'. (4.22) 



Proof. In view of Lemma 2.4 the claim readily follows from the general theory of 
exponentially low lying spectra for metastable Markov chains ||6|. For a direct proof 
note that under the assumptions of the Proposition, one (possibly after further 
shrinking the value of 5) can upgrade (4.20) as 



log max | (5^ + R X N 



XihNj(m) 
h Nj (mi) 



< -e, 



and (4.22) follows from the spectral theorem. 



(4.23) 
□ 



Asymptotics of Dirichlet Eigenvalues R^i- Define X(m) = \J\ — m 2 . The 

asymptotics of i% ; up to zero order terms is given [25J by 

= -JVri - Xo(m) + 



R\ 



N,l 



-Nri 



l V"(m{) 
X(m{) 







(4.24) 



where we used the explicit expression (4.4) for V in the second equality. Xo was 
defined in (4.16). The claim of Corollary 4.1 follows now from Proposition 4.2 □ 
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The Lagrangian Co was denned in (1.11) 

Lower bounds on Co- Fix a (z A and consider 0* = and, for j5 ^ a, Q\ — — ^ 
Recall that v G Rq, that is f a = — Ylp^a Vp. Therefore, for any a 



t. 



Co(m, v) > sup 



1 l fcOL 



m a mp \ a/ 3 (cosh(t) — 1] 



Define A Q (m) = Ylp^a yf m a m i3 For \v a \ > A Q (m) one may choose t* 

sign(f a ) log • We i n f er: If l^ol > A a (m), then 

\v n \ 



C Q (m,v) > \v a \ lo 



\ a (m) 



1 



(A.l) 



Upper bounds on the Lagrangian Co- Consider 

K-o{rn, v) = sup < ^ v a 9 a - ^ ^/m a m p A Qj(3 cosh (6^ - O 

Since C (m,v) = fco(m,v) + ^ Q A a (m), it would be enough to control the depen- 
dence of 7Z on v. 

Let us say that a flow / = {f a p} is compatible with v G Rq; / ~ t> if: 

(a) It is a flow: / a/3 = -fp a . 

(b) For any aei,^ fj3a = v a . 

Then £ v a 9 a = \ Y^aA ! 3 ~ <*)fo0- Hence, for any / ~ v, 



TZ = sup \ - ^(Op - 6 a )f a/3 - ^/m~m^ X a/3 cosh(0p - 6 a ) \ 



(A.2) 



We shall rely on the following upper bound on each term in (A.2): For any / and 
a > 

sup {ft -a cosh (t)} < |/| log ( 1 + — 
t \ a 

Consequently, we derive the following upper bound on TZ : 

I f a/3 1 



n (m,v)<mi VlMiog^i 



(A.3) 



Regularity of minimizers. Let m G int (A^). We claim that there exists 5q > 
and to > such that for any mf in the (^-neighbourhood of m the minimizer 7* of 

-to 

£(7(s),7(s))ds 



inf 

7(0)=m' ,y(t )=m Jo 

exists and is, actually, C°°. Indeed, an absolutely continuous minimizer exists by 
the classical Tonelli's theorem. By lower (A.l) and upper (A.3) bounds on the 
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Lagrangian, it is easy to understand that minimizers stay inside int(A^) once to an d 
Sq are chosen to be appropriately small. But then the regularity theory of either 11 
or |I] applies and yields Lipschitz regularity on [0,t ]. Since, the Lagrangian C is 
strictly convex in the second argument, and, in the interiour of A^, it is C°° in both 
arguments, the C°° of the minimizer follows from the implicit function theorem, see 
e.g. (7|. 
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